Frustrated magnets and quantum paramagnetic phases at finite temperature 
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We develop a general framework, which combines exact diagonalization in small clusters with 
a density matrix variational principle, to study frustrated magnets at finite temperature. This 
thermodynamic hierarchical mean-field technique is used to determine the phase diagram and mag- 
netization process of the three-dimensional spin-f/2 J1-J2 antiferromagnet on a stacked square 
lattice. Its non-magnetic phase exhibits a thermal crossover from a quantum to a classical para- 
magnet at a temperature T — To which can be extracted from thermodynamic measurements. At 
low temperature an applied magnetic field stabilizes, through order-by-disorder, a variety of phases 
with non-trivial spin textures and a magnetization plateau at half-saturation which continuously 
disappears at T ~ To. Our results are relevant for frustrated vanadium oxides. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.40.Cx, 75.60.Ej 



Introduction.- Frustrated magnetic materials have 
been focus of active condensed matter research in the 
past two decades^. In these systems, competing interac- 
tions and frustrated lattice topology often lead to fasci- 
nating effects, e.g. magnetic monopoles^^ in spin ice ma- 
terials or magnetization plateaux in an applied magnetic 
field^H and stabilize exotic quantum states of matter, 
such as spin liquids^ or valence bond solids^. Under- 
standing these paramagnetic phases, which do not break 
any obvious global symmetry and thus are not necessar- 
ily identified by an order parameter, is of fundamental 
interest in material physics. 

The theoretical description of frustration-driven phe- 
nomena is extremely challenging, especially in spatial di- 
mensions larger than one. For example, frustrating inter- 
actions render large-scale quantum Monte-Carlo (QMC) 
simulations impractical due to the "sign problem" 8 . 
Various approaches were developed to tackle frustrated 
magnets 9 . One class of methods proposes an expansion 
around a magnetically-ordered state (spin-wave and se- 
ries expansions), or in certain limiting cases, e.g. high 
temperature. Another class focuses on ground state 
(GS) properties of the system (coupled cluster and Lanc- 
zos methods). Meanwhile it is both theoretically and 
experimentally^ relevant to inquire: What are potential 
signatures of magnetic frustration in the thermodynamic 
properties of a given material? Clearly, quantum effects 
due to a non-trivial GS will become apparent at tem- 
peratures of the order of the characteristic energy scales 
involved in the formation of that particular GS. 

In the present Rapid Communication we address the 
above question. We develop an unbiased and gen- 
eral framework aimed at studying the interplay between 
quantum and thermal fluctuations in frustrated magnets. 
Our method couples the recently developed hierarchical 
mean-field (HMF) theory^ and the well-known thermo- 
dynamic variational principle 11 . A key idea is the real- 
ization that various competing local orders can only be 
captured within an exact diagonalization scheme, while 
transitions between them will be correctly described only 
in an infinite system. Hence, one starts by partitioning a 



lattice into relatively small spin clusters (degrees of free- 
dom) in accordance with point-group symmetries. These 
new degrees of freedom provide a language in which the 
original model Hamiltonian is represented. Approxima- 
tions are introduced in the form of a variational principle 
applied to the free energy in the new representation. Our 
approach relies on numerical as well as analytical efforts: 
provided the cluster is chosen properly, even a simple 
variational ansatz for the density matrix (DM) yields a 
complete phase diagram of the system. Results can be 
systematically improved by considering larger clusters or 
more complicated trial DMs. Once the system DM is 
known, any observable can be computed even inside non- 
magnetic (paramagnetic) phases, in contrast to the usual 
mean-field (MF) techniques which only yield instabilities 
of the magnetically-ordered states. 

We illustrate the thermodynamic HMF (THMF) for- 
malism by studying the phase diagram at finite temper- 
ature T and properties in an applied magnetic field of 
the "stacked" J1-J2 modet^tyj w hich describes a spin- 
1 /2 Heisenberg antiferromagnet on an orthorhombic lat- 
tice with first ( J\) and second ( J2) neighbor interactions 
in the a6-plane, and a nearest-neighbor (NN) exchange 
(Jc) along the c-axis. This model is a good approxi- 
mation for layered vanadium oxide material o 15 " 1 \ such 
as Li 2 VO(Si,Ge)OjiSland PbV0 3 19 . Previous works in- 
dicate that at T = the J\-J-2,-Jc model exhibits a 
quantum paramagnetic phase whose stability can be con- 
trolled by changing J c ■ We show that this phase persists 
even at finite T and introduce an important temperature 
scale, To, at which a crossover from quantum to classi- 
cal paramagnetic behavior takes place. Numerical value 
of T can be extracted from thermodynamic or magne- 
tization measurements. Below To, an interplay between 
quantum fluctuations inside the paramagnetic region and 
external field results in a variety of phases characterized 
by non-trivial magnetic orders, and leads to a magneti- 
zation plateau around half of the saturation field. Our 
findings are directly relevant for the frustrated perovskite 
PbV03 which shows no magnetic order and is believed 
to realize a J1-J2 quantum paramagnelP^l. 



The THMF method. - Let us consider a quantum spin 
model defined on a lattice with N sites and periodic 
boundary conditions. Following the HMF prescription^], 
we partition the lattice into Nu clusters of N q sites each, 
so that N — N\jN q . It is assumed that eigenstates of an 
isolated cluster are known. Each cluster state \a) can be 
associated with a Schwinger boson 7^, subject to a con- 
straint 7^7a = 1- It is important to note that the ver- 
sion of THMF method developed below treats this con- 
straint exactly, a condition that is crucial for the method 
to be variational. 

In terms of "cluster" degrees of freedom the original 
Hamiltonian of the system can be written as: 



H 



with Hq and V being the cluster self-energy and inter- 
cluster interaction, respectively. The subscript i labels 
different clusters in the coarse-grained lattice and sum- 
mation over doubly repeated indices is assumed. The 
range and type of couplings in the second term are de- 
termined by the original model. Each spin (and there- 
fore each cluster) has v\ links of type-A (A = 1, 2, . . .) 
with corresponding interactions V\. For instance, on a 
square J\-Ji lattice there are V\ = 2 NN and V2 = 2 
next-NN (NNN) links per site. The Hamiltonian 
operates on a Hilbert space spanned by the products 

|{°}) = rL7JajO)> where |0) is the unphysical Schwinger 
boson vacuum corresponding to all "empty" clusters. 

The THMF theory is a variational approach with re- 
spect to the free energy. Here we consider a simple trial 
DM Po [H MF ] = e~ K / T /Z with K = E^Hmf) ab jj a J lb , 
where the Boltzmann constant ks = 1, and the MF 
Hamiltonian Hmf, whose matrix elements play the role 
of variational parameters^, is self-consistently deter- 
mined by minimizing the free energy F = E — TS = 
Trp (H + T\np ) = Tr p (H - K) - T In Z . The parti- 
tion function can be expressed in terms of the eigenvalues 
En of H MF as Z = Z? D with Z x = J2n e ~ En/T - Then 
the free energy becomes: 

F 
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= — Tlog Zi+tiiUJi 



#0+E ^\^2(Vx)l2^2-H MF 



Here u = e ~ HMF/T ' /Z x is the single-cluster DM, "tr" 
denotes a trace over single-cluster MF states, and 

tr 1 tr 2 (Vx)i2WiW2 = {V\)l\llu aiali u a ^ 2 . 

A simple calculation yields the MF Hamiltonian^ 

(H MF )a'a = (H )a'a + 2 E ^(^^Mi ■ ( 2 ) 

A 

Using its eigenstates i?" one can compute the DM w: 

u ab =Y.*- En,TRn MT /E e_B " /T (3) 



and the free energy F/N = (l/2N q )Trcu(H - H MF ) - 
(T/N q )kiZi. In the limit T -> only the GS eigenpair 
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Figure 1. The J\-Ji-Jc lattice. Solid, dashed and dotted 
lines denote Ji , J2 and Jc exchange interactions respectively. 
Shaded cubes are clusters used in the THMF calculation. 



tr, 



(E , R°) contributes to the DM i.e. uj ab -> R°(R\ 
and we recover the T = HMF results^. 

Expressions ^ and ^ define the THMF self- 
consistent (in terms of Lo ab ) scheme. Any local observable 
(O) can be computed as (O) — trwC 

J\ -J2 -Jc model. - We now use the THMF method to 
study the J\-J%-Jc modeP 1 (Fig. [I): 

H = ( J i E + J 2 E + J c E ) S ^ - h zZ s ^ ( 4 ) 
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where all couplings Ji,i,c are positive, Si is a spin- 1/2 
operator at site i and h is an external magnetic field. This 
model describes a three-dimensional (3D) analog of the 
J\- J% antiferromagnelP^ defined on a cubic lattice with 
intralayer NN ((ij)) and NNN (((ij))) interactions J\ and 
J2 respectively, and an interlayer NN ((ij) z ) exchange 
Jc- In the following we adopt the units J\ = 1. 

At h = T = 0, the coupling Jc can be viewed as a 
tuning parameter controlling quantum effects associated 
with frustration. In particular, for Jc > Jc the quantum 
paramagnetic phase of the J1-J2 model disappears and 
the system exhibits a direct 1st order transition from 
Neel to a columnar antiferromagnetic (AF) state. For 
Jc < Jc there is an intermediate non-magnetic region 
which appears for a finite range of J 2 and vanishes at J®,. 
Thus the point [^(J^.), Jq] is a multicritical point where 
three phase boundaries converge. There seems to be a 
controversy regarding the order of phase transitions oc- 
curring at these boundaries: spin-wave studies^ predict 
one 2nd and two 1st order lines, while series expansions^ 
indicate that all phase boundaries are 1st order. Our 
analysis supports the spin- wave scenario^R 

We apply the THMF theory to compute the finite- 
temperature phase diagram of Hamiltonian Q. While 
the effect of thermal fluctuations on ordered phases is well 
known, their role inside a quantum non-magnetic state is 
unclear. Since this state does not break any continuous 
symmetry, the system can only undergo a crossover from 
a quantum to a standard classical paramagnet. We will 
show how the corresponding temperature scale can be 
deduced from experimentally accessible thermodynamic 
quantities, such as specific heat or uniform susceptibil- 
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Figure 2. Phase diagram of the J1-J2-JC model at h = 
(Ji = 1). There are three phases: (n, ir, 7r)-Neel antiferromag- 
net (AF), (0, n, 7r)-columnar AF (C), and paramagnet (P). 

ity. We also study properties of the J\-Ji-Jc model in 
an applied magnetic field. In the paramagnetic phase 
the magnetization curve has a plateau, similar to that 
of Ref. [55] obtained for the two-dimensional (2D) J1-J2 
model at T = 0. We demonstrate that with increasing 
field, a system exhibits a set of metamagnetic transitions 
characterized by non-coplanar spin textures. 

Application of the THMF technique starts by selecting 
a proper degree of freedom. Since the J1-J2-JC lattice 
has an orthorhombic unit cell with symmetry D2/J^, the 
smallest cluster compatible with this symmetry is a cube, 
shown in Fig. [I] Each cluster has 6 NN (4 in the afr-plane, 
2 along the c-axis) and 4 NNN clusters in the afr-plane; 
i.e. z/i = i>2 = 2 and vq = 1 in ([2|. Next, we compute the 
matrices H and V in Eqs. (JlJ and ^ using the method 
of Ref. [TUJ Finally, (jj a b is determined self-consistently 
using Eqs. ^ and (J3]> - 

Finite-T phase diagram.- The zero-field phase dia- 
gram of the J\-J2~Jc model Q is presented in Fig. [2] 
At T — and Jc < J c there exist three phases: AF with 
the wavevector Q = (tt, it, tt), plaquette quantum para- 
magnet (P), and columnar AF (C) with Q = (0, 7r, it) or 
(it, 0, 7r). The magnetic phases are characterized by an 
order parameter of the form Mi = (Si) = Mqc 1 ®™ 1 . The 
transition AF-P (P-C) is 2nd (1st) order for all values 
of Jc- This conclusion is consistent with the phase dia- 
gram of the J Y -J 2 (J c = 0) modePland agrees with the 
spin- wave^ and coupled cluster^ analysis. For Jc > Jc 
there is only a 1st order AF-C transition. The THMF 
calculation yields J^T = 0) ~ 0.28 - 0.30, in agreement 
with previous works^MLiJ. The real-space structure inside 
the P-region is plaquette crystal-like with layers covered 
with 2x2 plaqucttes, each being in its singlet GS. This 
state is similar to the paramagnetic GS of the 2D Ji- J 2 
modePSJ and is stabilized because for the P-phase Jc is 
quite small, so the system exhibits quasi-2D behavior. 

With increasing T the ordered AF- and C-states are 
suppressed via a 2nd order (classical) phase transition 
at a Neel temperature Tjy(J2, Jc)i which is accompa- 
nied by a jump in the magnetic specific heat C v , see 
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Figure 3. (a) Specific heat C v for the phases in Fig. [2] at 
Jc = 0.1: AF (J 2 = 0.1), P (J 2 = 0.5) and C (J 2 = 0.8). (b) 
Uniform susceptibility xo for the paramagnetic phase (main 
plot), and magnetic states (inset) with arrows indicating the 
Neel temperature. Parameters are the same as in panel (a). 

Fig. [3ja). Since the THMF method ignores long-range 
fluctuations, Tjy remains finite even in the 2D limit 
Jc — 0, which constitutes a violation of the Mermin- 
Wagner theorerrPS common to MF theories. We note 
that in the C-state along with the columnar magnetiza- 
tion M c one can introduce an Ising order paramete r^ 5 * 26 ! 

(7 = (l/N)Y^ Xi y S Xt y(S X+ l >y - S x> y + l), Wltll tllC SUHl- 

mation extending over N sites of the original lattice and 
(x,y) = i. This Z2 order parameter is not subject to the 
conditions of the Mermin- Wagner theorem and vanishes 
at a non-zero temperature 2 ^ via a 2nd order phase tran- 
sition. However, at the HMF level a and M c vanish at 
the same Neel temperature. 

Inside the P-phase, a crossover takes place from quan- 
tum (at low T) to classical (at high T) behavior. Because 
the paramagnetic GS is gapped, C v exhibits a peak [Fig. 
|3][a)] . We argue that the position of this peak serves as an 
estimate of the crossover temperature To at which ther- 
mal fluctuations become comparable to the gap in the 
GS. For instance, at J 2 = 0.5 and J c = 0.1, T - 0.3. 
The gap in the P-state manifests itself in an activated 
behavior of the uniform linear magnetic susceptibility 
Xo(T) = dM z /dh\h=o, presented in Fig. [3jb) . Contrary 
to the ordered AF- and C-states where Xo(T — 0) is fi- 
nite and has a peak approximately at the corresponding 
Tjv, for the P-phase Xo{T — 0) = and shows a broad 
maximum around To- Since the P-phase has a plaquette 
structure, it is natural to examine the behavior of various 
plaquette "order parameters". We considered the func- 
tions f 4 - (i/A0 J2 x , y Sx, v [(-iys x +i, y + {-i)ys Xt y +l \ 

and Q = (1/2N) £(^1234 + ^1234) introduced in Refs. M 
and 1271 respectively. In the expression for Q the summa- 
tion takes place over plaquettes in a&-planes and P1234 is 
an operator of cyclic permutation of plaquette vertices. 
Both functions remain finite and decay as ~ 1/T at large 
temperature. The crossover manifests itself through a 
peak in dF^/dT and dQ/dT around Tq. As any real-space 
method, the THMF theory involves explicit translational 
symmetry breaking (cf. Ref. I2ip . predicting a crossover 
inside the paramagnetic phase. In an exact thermody- 
namic limit solution this crossover may become a phase 
transition because of melting of the plaquette crystal. 

Now let us consider transitions between different 
phases, which are triggered by tuning J2 while keeping 
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Figure 4. (a) Magnetization process of the J\-Ji-Jc model 
with J 2 = 0.5, Jc = 0.1, and M aat = 1/2. Main panel: Total 
magnetization M versus applied field at T = and schematic 
spin profiles in the corresponding field ranges (the length of 
the arrows is proportional to the magnitude of the spin ex- 
pectation values). Inset: M at finite temperature T = 0.16 
(solid line) and T — 0.3 (dotted line), (b) and (c) Tempera- 
ture dependence of the specific heat and M at ft = 2.5. The 
2nd order phase transition happens at T ~ To. 

T > and Jc fixed. The transition AF-P remains 2nd 
order. On the other hand, the transition P-C is clearly 
discontinuous at low T because symmetries of the two 
phases are not related by the group-subgroup relatiorP^l 
At higher T (still not destroying the magnetic order), 
the jump in the columnar magnetization vanishes and 
the transition becomes continuous. 

Finally we compare the critical temperature and expo- 
nents obtained from the 8-spin cluster THMF with their 
known values. For the unfrustrated 3D Heisenberg model 
(J 2 = and Jc = 1) the specific heat and order param- 
eter exponents are o>thmf — and Pthmf — 0.470(1), 
and the Neel temperature is T^ HMF = 1.308. These val- 
ues should be compared with the results of the Weiss 
molecular field aw = 0, /3w — 0.5, = 1.5, and 
QMCPHH p QMC = o.36, T® MC = 0.946. Near the para- 
magnetic phase, e.g. for J 2 = 0.3 and Jc =0.1, we have 
Pthmf = 0.474(2) and T FHMF = 0.512. The THMF 
exponents and Tjy can be improved by increasing the 
cluster size and implementing a finite-size extrapolation. 

Magnetic states in an applied field.- An exhaustive 
spin-wave study of the magnetization process M(h) in- 
side ordered phases of the J\-Ji-Jc model was performed 



in Ref. [THl On the contrary, high-field properties of the 
quantum paramagnetic state received much less atten- 
tion, with efforts exclusively focused on the case Jc — 
(2D Jx-Ji model). Particularly, in Ref. [22] a half- 
saturation [M = 1/4) magnetization plateau character- 
ized by a collinear spin ordering was proposed around the 
maximally frustrated point J 2 = 0.5. 

Here we apply the THMF theory to study field-induced 
metamagnetic transitions inside the non-magnetic region 
(Fig. [2]) of the J\-J%-Jc model with J 2 = 0.5 and 
Jc = 0.1. The main panel of Fig. |4]ja) displays the T = 
magnetization curve. For h below certain threshold value 
fiMi the magnetization vanishes due to the spectral gap. 
For h > h,M, the system exhibits a set of phases with non- 
coplanar magnetic textures, shown in the figure. While 
the spin orderings at small and large fields simply reflect 
canted AF sublattices, the structures immediately before 
and after the 1/2-plateau are non-trivial because the clas- 
sical canting angle varies for different spins. The plateau 
state has long-range order characterized by a collinear 
magnetic structure with two spins per cluster antiparal- 
lel to the field (cf. Ref. l22j) and an Ising order parameter. 

The plateau width [Fig. |4][a), inset] and threshold field 
Km vanish at T ~ To in agreement with the behavior of 
C v [Fig. [3]ja)]. For a fixed field inside the plateau, the 
magnetic order collapses via a 2nd order phase transition 
at T ~ T . In Fig. ^b) and (c) we show the temperature 
dependence of C v and M at h = 2.5. The non-monotonic 
behavior of M(T) is easy to understand by observing 
that at T = spins are "locked" in a specific pattern by 
the interactions. Thermal fluctuations unlock the spins 
allowing them to orient along the field, thus increasing 
M(T) before the transition. This highly non-trivial mag- 
netization process is specific to the P-state in Fig. [2j and 
we propose to use it in conjunction with the peak in C' v 
as characteristic signatures of a quantum paramagnet. 

Conclusion.- We developed and applied the THMF 
method to address the interplay between quantum and 
thermal fluctuations in the spin-1/2 Jx-J-x-Jc antiferro- 
magnet. Focusing on the non-magnetic region of the 
model, which is inaccessible for other theoretical tech- 
niques, we studied the crossover between a classical (due 
to thermal effects) and a quantum paramagnet, and 
demonstrated how the crossover temperature scale To can 
be extracted from thermodynamic and high-field mea- 
surements. At low temperature T < T) quantum fluc- 
tuations inside the paramagnetic state are manifested in 
a variety of field-induced spin structures and a magne- 
tization plateau at half-saturation. Our results can be 
verified in experiments with vanadium oxides of the type 
PbV0 3 . AssumingSH that J x ~ 70 - 100 K (and J 2 , J c 
inside the paramagnetic phase), one gets To ~ 20 — 30 K 
and hu ^ 25 — 35 T. The magnetization plateau should 
become apparent for h ~ 100 — 200 T. 
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